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o ■ 

X) ■ ABSTRACT 

l/^ ' We report cosmic-ray proton and helium spectra in energy ranges of 1 to 120 GeV 

^ ■ and 1 to 54 GeV/nucleon, respectively, measured by a balloon flight of the BESS spec- 

trometer in 1998. The magnetic-rigidity of the cosmic-rays was reliably determined by 
highly precise measurement of the circular track in a uniform solenoidal magnetic field 
I of 1 Tesla. Those spectra were determined within overall uncertainties of ±5 % for 

I protons and ±10 % for helium nuclei including statistical and systematic errors. 

o 
o 
o 

o 

^ ' The absolute fluxes and spectra of primary cosmic-ray protons and helium nuclei are funda- 

mental information as references in cosmic-ray physics. They are needed to calculate secondary 
anti-protons, positrons, and diffuse gamma radiation, which in turn provide important knowledge 
of particle propagation and the matter distribution in interstellar space. Those are also indispens- 
5^ i able for studying atmospheric neutrinos. Although measurement of the proton and helium energy 

spectra has been performed in various experiments, their resultant absolute fluxes show discrepan- 
cies up to a factor of two at 50 GeV. This ambiguity causes large uncertainty in calculations of the 
atmospheric neutrinos, as well as secondary anti-protons, positrons, and diffuse gamma-rays. 
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We report a new precise measurement of the cosmic-ray proton and helium spectra over the 
energy ranges of 1 to 120 GeV for protons and 1 to 54 GeV/nuclcon for helium nuclei, based on a 
half of the data from a BESS balloon flight in 1998. The covered energy range is relevant to the 
atmospheric neutrinos observed as "fully contained events" in Super-Kamiokande. In the BESS-98 
flight, a new trigger mode was implemented with a silica-aerogel Cherenkov counter to record all 
energetic particles instead of sampling the protons at a ratio of 1/60 as done in the previous BESS 
flights. This drastically improved statistics in the high-energy region above 6 GeV/nucleon, as 
reported here. 

2. THE BESS SPECTROMETER 

The BESS detector is a high-resolution spectrometer with a large acceptance to perform highly 
sensitive searches for rare cosmic-ray components, as well as precise measurement of the absolute 
fluxes of various cosmic-ray particles (Orito 1987; Yamamoto et al. 1994; Ajima et al. 2000). As 
shown in Figure 1, all detector components are arranged in a simple cylindrical configuration with 
a thin superconducting solenoidal magnet. In the central region, the solenoid provides a uniform 
magnetic field of 1 Tesla (±7 % in a fiducial volume). The trajectory of an incoming charged particle 
is measured by using the tracking system which consists of a central JET chamber and two inner 
drift chambers (IDC's) in a volume of 0.84 ni in diameter and 1 m in length. The magnetic-rigidity 
{R = pc/Ze) is reliably determined by a simple circular-fitting of the deflection (i?^^) using up-to 
28 hit points each with a spatial resolution of 200 fim . Figure 2 shows the deflection uncertainty 
for protons evaluated in the track fltting procedure. Since the magnetic field is highly uniform, 
it has a narrow and sharp peak around 0.005, which corresponds to maximum detectable rigidity 
(MDR) of 200 GV, and has a small tail. The outermost detector is a set of time-of-flight (TOE) 
hodoscopes with 2 cm thick plastic scintillators. It provides the velocity {(5 = v/c) and energy loss 
{dE/dx) measurements. The time resolution for energetic protons in each counter is 55 ps rms, 
resulting in a resolution of 1.4 %. The data acquisition sequence is initiated by a first- level TOF 
trigger by a simple coincidence of signals in the top and bottom scintillators, with the threshold 
level of 1/3 pulse height from minimum ionizing particles (MIP's). If the pulse height exceeds 2.5 
times MIP signal, the TOE trigger is labeled as "helium-trigger," otherwise "proton-trigger." The 
TOF trigger efficiency was evaluated to be > 99.95 % from the muon data taken at sea level; thus, 
the systematic error caused by the TOF trigger inefficiency is negligibly small. In order to build 
a sample of unbiased triggers, one of every 60 and 25 "proton-triggered" and "helium-triggered" 
events, respectively, was recorded irrespective of succeeding online-selections. 

The instrument has a threshold-type Cherenkov counter with a silica-aerogel radiator just 
below the top TOF hodoscope (Asaoka et al. 1998). The radiator was newly developed prior to the 
BESS-98 flight and it has a refractive index of 1.022. An auxiliary trigger was generated by a signal 
from the Cherenkov counter to record all of the high energy particles above ~ 6 GcV without bias 
nor sampling. The efficiency of the Cherenkov trigger for the energetic particles was determined to 
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be 92.1 lb 3.0 % by using the data sample of unbiased triggers. 

The BESS spectrometer was flown 1998 July 29-30 from Lynn Lake, Manitoba, Canada. It 
floated at an altitude of 37 km (residual atmosphere of 5 g/cm^) with a cutoff rigidity of 0.5 GV 
or smaller. The solar activity was close to the minimum. 

3. DATA ANALYSIS 

In the first stage of data reduction, we selected events with a single track fully contained inside 
a fiducial volume defined by the central four columns out of eight columns in the JET chamber. This 
definition of the fiducial volume reduced the effective geometrical acceptance down to ~ 1/3 of the 
full acceptance, but it ensured the longest track-fitting and thus the highest resolution in the rigidity 
measurement. The single-track selection eliminated rare interacting events. In order to verify the 
selection, events were scanned randomly in the unbiased trigger sample, and it was confirmed that 
995 out of 1,000 visually-identified single-track events passed this selection criteria and interacting 
events were fully eliminated. Thus, the track reconstruction efficiency was 99.5 =b 0.2 % for a 
single-track event. In order to assure accuracy of the rigidity measurements, event quality such 
as in the track fitting procedure was imposed on the single-track events. The efficiency of this 
quality-check was estimated from loose cuts of the flight data to be 93.8 =b 0.3 % and 93.0 it 0.9 % 
for protons and helium nuclei, respectively. It was almost constant over the whole energy range. 

In the final stage of data reduction, particle identification was performed as shown in Figure 3 
by requiring proper dE/dx and 1//3 as a function of rigidity. According to a study of another 
sample of 5 x 10^ protons and 4 x 10^ helium nuclei selected by using independent information of 
energy loss inside the JET chamber, the dE/dx selection efficiencies were 99.3 ± 0.2 % for protons 
and 98.2 ±0.7 % for helium nuclei, and the contamination probabilities should be less than 3 x 10~^ 
for protons and 4 x 10~^ for helium nuclei. Since the 1/(3 distribution is well described by Gaussian 
and the half width of 1//3 selection band was set at 3.89 cr, the efficiency is very close to unity 
(99.99 % for pure Gaussian). 

Very clean proton samples were obtained below 3 GV. However, deuterons started to contam- 
inate the proton band around at 3 GV, where the contamination was observed to be 2 %. No 
subtraction was made for this contamination, because it was as small as the statistical errors and 
a deuteron-to-proton ratio decreases with energy (Sco ct al. 1997) following a decrease in escape 
path lengths of primary cosmic-ray nuclei (Engelmann et al. 1990). In conformity with previous 
experiments, all doubly charged particle were treated as ^He. With the data reduction described 
above, 826,703 protons and 77,325 helium nuclei were finally identified. The combined efficiencies 
were 92.7 it 1.0 % for protons and 90.8 it 1.5 % for helium nuclei. 

The geometrical acceptance defined for this analysis was calculated to be 0.0851 it 0.0003 m^sr 
for energetic particles by using simulation technique (Sullivan 1971). The simple cylindrical shape 
and the uniform magnetic field make it trivial to determine the acceptance precisely. The error 
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axises from uncertainty of the detector alignment within 1 mm. The ratio of hve data-taking-time 
was measured exactly to be 86.4 % by counting 1 MHz clock pulses. 

The energy of each particle at the top of the atmosphere can be calculated by summing up 
the ionization energy losses with tracing back the event trajectory. Detection efficiencies were 
studied by using Monte Carlo (M.C.) simulation. The M.C. code was developed to incorporate 
detailed description of various interactions of helium nuclei into GEANT (Brun et al. 1994), where 
the cross-sections and angular-distributions of the nuclear interactions were evaluated by fitting 
experimental data (Bizard et al. 1977; Abdurakhimov et al. 1981; Gasparyan et al. 1982; Ableev 
et al. 1985; Grebenjuk et al. 1989; Glagolev et al. 1993; Abdullin et al. 1994) to energy dependent 
empirical formulas (Bradt & Peters 1950). The electromagnetic processes, mainly due to (5-rays, 
are also treated properly. They are more significant in helium nuclei interactions, because cross- 
sections in electromagnetic processes (a"em) behaves as ~ whereas those in hadronic processes 
(o'had) is approximately proportional to (2Z)^/^. The M.C. well reproduced the observed event 
shape. The simulated and observed number of hit-counters in the bottom TOP hodoscope, for 
instance, were agree within a discrepancy of 0.9 % and 1.8 % for proton and helium, respectively. 
The systematic errors in the M.C. originate mainly in uncertainties of (Jhad and (Jem- We attributed 
relative errors of ±5 % to (Jhadip + A), ±5 % to (Tem(p + A), ±15 % to o-/iad(He + A), ±20 % to 
c7em(He ± A), and ±20 % to (j/iad(CNO ± A). Another source of a systematic error in the M.C. was 
the uncertainty of the material distribution inside the BESS spectrometer, which was estimated 
to be ±10 %. The probability that cosmic-ray protons and helium nuclei, respectively, can pass 
through the whole detector without interaction was 87.6 ± 2.3 % and 74.5 ± 6.8 % at 1 GeV, and 
79.7 ± 2.9 % and 64.6 ± 7.5 % at 100 GeV. 

According to similar M.C. studies, the probabilities that primary cosmic-rays can penetrate the 
residual atmosphere of 5 g/cm^ is about 94 % and 90 % for the proton and helium, respectively, over 
the entire energy range discussed here. Atmospheric secondary protons, which account for about 
3.5 % at 1 GeV and less than 1.5 % above 10 GeV of observed protons, were subtracted based 
on the calculation by Papini et al. (1996). Atmospheric secondary helium above 1 GeV/nucleon is 
dominated by fragments of heavier cosmic-ray nuclei (mainly Carbon and Oxygen) . The flux ratio 
of atmospheric secondary helium to primary C+0 was calculated to be 0.14 at a depth of 5 g/cm^, 
based on the total inelastic cross-sections of CNO ± Air interactions and the helium multiplicity 
in ""^^C + CNO interactions (Ahmad et al. 1989). The total correction of atmospheric secondary 
helium due to all nuclei with Z > 2 was estimated to be about 2 % over the entire energy range. 
The ambiguity in this estimation arises from the uncertainties in ahadi^^O + A) and in absolute 
fluxes of heavier cosmic-ray nuclei, to which we attributed relative error of ±20 %. The systematic 
errors originating in the correction of the residual air effect were estimated to be ±0.3 % for protons 
and ±2.0 % helium nuclei. 

Because of the finite resolution in rigidity measurement, and the very steep spectral shape, 
the observed spectrum shape may suff'er deformation. The effect of finite resolution was estimated 
by simulation, where the error in rigidity measurement was tuned to reproduce the distribution 
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shown in Figure 2. The effect was found to be smaller than 1 % below 25 GV, but it became visible 
with increasing rigidity. The observed spectrum gradually gets lower than original spectrum, with 
a ratio of —2.5 % at 70 GV and then rapidly rises to ±0 % around 120 GV. No correction was 
made for this deformation, because the effect is as small as the statistical errors over the energy 
range discussed here. 



4. EXPERIMENTAL RESULTS AND CONCLUSION 

The proton and helium fluxes at the top of the atmosphere have been obtained from the BESS- 
98 flight data as summarized in Table 1 and as shown in Figure 4 in comparison with the previous 
experiments. The first and second errors in Table 1 represent statistical and systematic errors, 
respectively. The overall errors including both errors are less than ±5 % for protons and ±10 % 
for helium nuclei. The dotted lines in Figure 4 indicate the spectra assumed in the calculation of 
atmospheric neutrinos (Honda et al. 1995). Our results, as well as other recent measurements, are 
favorable to lower fluxes than the ones assumed in the atmospheric neutrino calculation especially 
above a few tens of GeV. It may suggest an importance of the reconsideration for the atmospheric 
neutrino flux predictions. Precise measurements of primary cosmic-rays will help to improve the 
accuracy in the atmospheric neutrino calculations. 

The authors would thank NASA and NSBF for the balloon flight operation. This experiment 
was supported by Grants-in-Aid from Monbusho and Hciwa Nakajima Foundation in Japan and by 
NASA in the U.S.A. The analysis was performed with the computing facilities at ICEPP, University 
of Tokyo. 
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Table 1. Proton and helium fluxes at the top of the atmosphere. 



Proton Helium 



Energy Range Flux ± AFluxgta ± AFluxgyg Ek Flux ± AFluxsta ± AFluxsy 

[GeV/n] [GcV/n] [ur'^sr-'^scc-^GcV /n)''^] [GeV/n] [rn-2sr-^sec-^(GoV/n)-^] 



1.00 


1.17 


1.08 


8.92 


± 


0.12 


± 


0.22 


xlO^ 


1.08 


8.21 ± 0.26 ± 


0.62 


xlO 




1.17 


1.36 


1.26 


7.72 


± 


0.11 


± 


0.19 


xlO^ 


1.26 


6.57 ± 0.22 ± 


0.50 


xlO 




1.36 


1.58 


1.47 


6.74 


± 


0.09 


± 


0.17 


xlO^ 


1.47 


5.46 ± 0.19 ± 


0.41 


XlO 




1.58 


1.85 


1.71 


5.46 


± 


0.08 


± 


0.14 


xlO^ 


1.71 


4.38 ± 0.15 ± 


0.33 


XlO 




1.85 


2.15 


2.00 


4.52 


± 


0.07 


± 


0.11 


xlO^ 


2.00 


3.29 ± 0.12 ± 


0.25 


XlO 




2.15 


2.51 


2.33 


3.63 


± 


0.05 


± 


0.09 


xlO^ 


2.33 


2.69 ± 0.10 ± 


0.21 


XlO 




2.51 


2.93 


2.71 


2.83 


± 


0.04 


± 


0.07 


xlO^ 


2.71 


1.90 ± 0.08 ± 


0.15 


XlO 




2.93 


3.41 


3.16 


2.22 


± 


0.04 


± 


0.06 


xlO^ 


3.15 


1.38 ± 0.07 ± 


0.11 


XlO 




3.41 


3.98 


3.68 


1.71 


± 


0.03 


± 


0.05 


xlO^ 


3.68 


1.12 ± 0.05 ± 


0.09 


XlO 




3.98 


4.64 


4.30 


1.27 


± 


0.02 


± 


0.03 


xlO^ 


4.30 


8.65 ± 0.44 ± 


0.66 






4.64 


5.41 


5.01 


9.65 


± 


0.19 


± 


0.26 


xlO 


5.01 


5.80 ± 0.34 ± 


0.44 






5.41 


6.31 


5.84 


6.89 


± 


0.15 


± 


0.19 


xlO 


5.84 


4.27 ± 0.27 ± 


0.33 






6.31 


7.36 


6.81 


4.91 


± 


0.02 


± 


0.20 


xlO 


6.80 


2.96 ± 0.04 ± 


0.24 






7.36 


8.58 


7.94 


3.43 


± 


0.01 


± 


0.14 


XlO 


7.92 


1.99 ± 0.03 ± 


0.16 






8.58 


10.0 


9.25 


2.42 


± 


0.01 


± 


0.10 


XlO 


9.24 


1.44 ± 0.03 ± 


0.12 






10.0 


11.7 


10.8 


1.70 


± 


0.01 


± 


0.07 


XlO 


10.8 


9.98 ± 0.20 ± 


0.81 


XlO" 


1 


11.7 


13.6 


12.6 


1.18 


± 


0.01 


± 


0.05 


XlO 


12.6 


6.82 ± 0.15 ± 


0.55 


xlO" 


1 


13.6 


15.8 


14.6 


8.05 


± 


0.04 


± 


0.33 




14.6 


4.50 ± 0.12 ± 


0.36 


xlO" 


1 


15.8 


18.5 


17.1 


5.57 


± 


0.03 


± 


0.23 




17.1 


3.16 ± 0.09 ± 


0.26 


XlO" 


1 


18.5 


21.5 


19.9 


3.78 


± 


0.03 


± 


0.16 




19.9 


1.99 ± 0.07 ± 


0.16 


XlO" 


1 


21.5 


25.1 


23.2 


2.51 


± 


0.02 


± 


0.10 




23.2 


1.51 ± 0.05 ± 


0.12 


xlO" 


1 


25.1 


29.3 


27.1 


1.67 


± 


0.01 


± 


0.07 




27.0 


9.15 ± 0.39 ± 


0.74 


XlO" 


2 


29.3 


34.1 


31.5 


1.10 


± 


0.01 


± 


0.05 




31.5 


5.98 ± 0.29 ± 


0.49 


XlO" 


2 


34.1 


39.8 


36.8 


7.35 


± 


0.08 


± 


0.31 


xlO"^ 


36.9 


4.30 ± 0.23 ± 


0.35 


xlO" 


2 


39.8 


46.4 


42.9 


4.87 


± 


0.06 


± 


0.20 


xlO"^ 


42.9 


2.65 ± 0.17 ± 


0.22 


XlO" 


2 


46.4 


54.1 


50.0 


3.22 


± 


0.05 


± 


0.14 


xlO"^ 


49.8 


1.88 ± 0.13 ± 


0.16 


XlO" 


2 


54.1 


63.1 


58.3 


2.10 


± 


0.04 


± 


0.09 


xlO"^ 












63.1 


73.6 


67.9 


1.36 


± 


0.03 


± 


0.06 


xlO"^ 












73.6 


85.8 


79.3 


9.17 


± 


0.20 


± 


0.39 


xlO"^ 












85.8 


100. 


92.6 


6.08 


± 


0.15 


± 


0.26 


xlO-2 












100. 


117. 


108. 


4.00 


± 


0.12 


± 


0.17 


xlO"^ 
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Fig. 1. — Cross-sectional view of the BESS instrument. 
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Deflection Uncertainty [GV" ] 



Fig. 2. — Deflection uncertainty for protons. Each area of the histogram is normahzed to unity. 
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Fig. 3. — Proton bands in (a) dE/dx (top-TOF) vs rigidity plane; and (b) 1/(5 vs rigidity plane 
after proton dE/dx selection. dE/dx in the bottom-TOF is also checked. The superimposed graphs 
show the proton selection criteria above 10 GV. Helium nuclei were selected in the same manner. 
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Fig. 4. — Absolute differential proton and helium spectra. Filled squares show results of the 
BESS-98 experiment. The spectra obtained by other experiments are also shown by different 
symbols indicated in the figure. Dashed lines show assumed spectra in the atmospheric neutrino 
flux calculation. 



